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^ ' In this work, we continue the study of the supersymmetric t-J model with 1/r^ 

(»— ^ ' hopping and exchange without translational invariance. A set of Jastrow eigen- 

^^ ' functions are obtained for the system, with eigenenergies explicitly calculated. The 



ground state of the t-J model is included in this set of wavefunctions. The spectrum 
of this t-J model consists of equal-distant energy levels which are highly degenerate. 
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In recent years, there have been considerable interests in study of low dimensional elec- 
tronic models of strong correlation, due to the possibility that the normal state of the two 
dimensional novel superconductivity [|I| may share some interesting feature of ID interact- 
ing electron system (non- Fermi liquid behavior). The one-band two dimensional Hubbard 
model reduces to the t-J model in the large on-site energy limit. The Hubbard model and 
the t-J model have been under intense study through various approaches. For these strongly 
correlated electron models, few exact results may be obtained in two dimensions. The high 
dimensional versions are much harder to study than their one dimensional ones. In one 
dimension, however, the Bethe-ansatz technique may allow us to exactly solve Hamiltonians 
in some special cases, such as the Lieb-Wu solution and the ordinary t-J model at its 
supersymmetric point 0. In particular, the ID long range exactly solvable electronic mod- 
els have attracted a lot of attention, since they display interesting physics with solutions of 
simple mathematical structure |J|,|,0,|,|,|ig,|n],|g,|l|,|l|,|15|,|l|,|l3,|l|,|l|,|O|,|ll. 

Recently, we have been able to explicitly construct all the constants of motion for the 
translational invariant long range supersymmetric t-J model, by mapping the system to a 
mixture of fermions and bosons, with the super-algebra representation for the electron fields 
| |12| |. Moreover, we have introduced a one dimensional supersymmetric t-J model with 1/r^ 
hopping and exchange without translational invariance. This system has also been shown 
by us to be completely integrable, with infinite number of conserved quantities explicitly 
constructed |jl2[. In this work, we continue the study of this t-J model. A set of Jastrow 



eigenfunctions, as well as their eigenenergies, are obtained explicitly. The ground state of 
the t-J model is included in this set of wavef unctions. We also briefly discuss the structure 
of the full spectrum for the system. 

The Hamiltonian for the one-dimensional t-J model is given by ||12| 

N 



Htj = {1/2)Pg 






P< 



G, 



where we take the hopping matrix and the spin exchange interaction to be tij/2 = Jij = 
l/{ri — rjY . L is the number of sites on the chain. In this model, the positions of the sites 



{vi} are given by the roots of the Hermite polynomial Hl{x) , and the spin component a takes 
values from 1 to A^. The operator q^ is the the creation operator for an electron at site i with 
spin a, Cicr is the corresponding annihilation operator. Their anti-commutation relations are 
given by {q^^, c]^J+ = 5,^5^,^^., {q^,, CjvJ+ = 0, {c^, c]^J+ = 0. The Gutzwiller projection 
operator Pq projects out all the double or multiple occupancies, Pq = Ili=i Poi'i^) , and 

-Pg(^) = So,n, + (^l,n,, with rii = J2a=l4aCia- The operator Pij = J2a=lT.a' = l4aCia'c]^,Cja 

exchanges the spins of the electrons at sites i and j, if both sites are occupied, and is zero 
otherwise. At half-filling, our t-J model becomes the long range spin model, introduced first 



by Polychronakos on such a non-translational- invariant lattice |14 



In terms of the b and / fields, the eigenequation of the t-J model can be written as [12 



- E {q^-q,)-'M,,(|,{{q},{a})=E<P{{q},{a}), (2) 

where ^({g}, {a}) = 0(gi<Ti, 520-2, ■ ■ ■ , qN,o-NjqNe+i, (1n,+2, ■■■ .qi) is the amplitude for the 
/ fermions to be at gi, q2,- ■ ■ , qwe^ while the spinless b bosons are at qNe+i^ qNe+2, ■ ■ ■ ,qL- 
Here, {a} = (ai, (T2, ■ ■ ■ , (XivJ and {q} = (gi, g2, ■ " " , ^l) = {xi,X2,- ■ ■ ,XN,,yi,y2r ■ ■ ^Vq)- 
The wavefunction is symmetric in the b boson positions {y}, while antisymmetric when 
exchanging Xiai and XjCTj. The operator Mij exchanges only the position variables qi and 
qj, that is, Mij(f){{q},{a}) = 0({g'},{a}), with {q} = (gi, ^2, ■■■,?*,■■■, ^j, ■■■, ^l) and 
{q'} = (gi) q2,- ■ ■ ,qj,- ■ ■ ,qi,- ■ ■ , qb)- in this approach, the / fermions and the b bosons 
occupy the whole chain, i.e. {g} and {g'} are permutations of the sites {ri, r2, ■ ■ ■ , r^}, and 
we work in the Hilbert space where at each site there is exactly one particle. Q, the number 
of the b bosons, is also the number of holes in the original problem; iVg, the number of the 
/ fermions, is also the number of the c electrons on the lattice. A^^, with a = 1,2,- ■ -N, 
the number of the / fermions with spin component a, is also the number of the c electrons 
with spin component a on the chain. 

The Hamiltonian in the first quantization, as given by Eq. (^, is 

H = - E iq.-q,r'M,^. (3) 



It commutes with the permutation operator Tij = P^jMij exchanging the / fermion spin and 
position simultaneously. Let us work in the Hilbert space where the number of fermions of 
each flavor is fixed, i.e. N^r, a = 1,2, ... ,N, is fixed. Consider the following wavefunction 
in Jastrow product form, 



[XiCri,X2<T2,- ■■ ,XN,CrNjyi,y2, 
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— T .\^''t''3P^^^9n(cT^~Uj) 



Jb f iXj -J 



j'e 2 



(4) 



where {x} and {y} span the whole lattice. We would like to show that this wavefunction is 
an eigenstate of the system. The Hamiltonian in Eq. (0) can be broken up into three parts, 
the first part Hi exchanges the / fermions, the second H2 exchanges the h bosons, and the 
third H^ exchanges the bosons and the fermions: 

^i = (-i) E {q^~q,VM,^ 

l<i<j<Ne 

H2 = (-1) E (?" - ?/3)"'M„/3 

l+Afe<a</3<i 

^3 = (-i) E E (g«-g.)-'M„,- (5) 

Afe+l<a<L l<i<A^e 

We then calculate the effects of these three parts when acting the Jastrow wavefunction 
given by Eq. (HI). The contribution from H2 is immediate: 



H2(t) = - E (?/° - yp) 

a</3 



(6) 



The contributions from Hi and H-^ are harder to deal with since many particle terms are 



involved. Using a similar trick introduced in Refs. 18,12], we have 



where 



^3</>=-EE(^^-2/")"'n 



,ya ^j\5a 



i a 



ji^i) ^'' ^^ 



-EE(^^-?/-) ' 11(1 + ^-.- 



ya -^i 



yya -^ijy-^i •^j) 



— Rest, 



N N, 



Rest = E E E E [J2(y» - ^* 



,r-3 



n 



a a=l r=3 JCp'^ jgj 

|J|=r 



XjeXj/xi 



yXi '^jj 



(7) 



(8) 



with p"" = {A; G {1, 2, ■ ■ ■ , A^e}; Cfc = cr}, Xj = {xj;j G J}. Then using the fact that for any- 
set X = (xi, 0:2, ■■■ , Xn), we have the following identity ( see Appendix A ) 

1 






0, 



(9) 



for all t = 0, 1, 2, ■ ■ ■ , ra — 2, we conclude that Rest = 0. 

The contribution from Hi is calculated in a similar manner: 



i^i0 = -E 



Xj Xj) 



;i-2<5...,) n 



,XA: •'^j \5a.a, /-^k ^i\S, 



■j'^k . 



i<j Ht^iJ) ^^ ^* ^'^ ^^ 



5I^^,<x,(a;i-Xj) 20-E(xi-Xj) 2(i_(5^^^j J| (1 + 5, 



£ — A -L 



Kj 



t<J 



Ki^hj) 



Xj^ Xi 



Xj^ Xj 



Y. '^^«<^. i^i - ^j) ^(P -Y(^i- 






i¥=j 



t<j 



Xi X^ } 11 Orr^rr^ I 



Kj 



E n^^. 






XiC{l,2,--,Afe}/ii k&Kx 

K2(l{l,2,-,Ne}/ij 

KiUK2^0 



Xk Xi 



n 

k€K2 



Xf^ Xj 



(10) 



Using the sum rule Eq. (H), and the fact that J^i^tj^tki^k — Xi) ^{xk — Xj) ^Sa^^jSaiUk = 0, 
the last term in the above equation becomes 



-0E(l-'^a,;<7,) E [i^i-Xj) ^{xk-Xi) ^6^^^^ - {xi - Xj) ^{xk-Xj) ^6^.^^ 

i<J k(¥=i,j) 

= -^Y. E i.Xi-Xj)~^{,Xk-Xi)-^5^^^^. 
i^j k(^i,j) 

In the end, we have 



(11) 



H(l) 



E (^^ 

i<i<j<L 



Qj^ 



+ E(^' - 


-Xj) ^6^^^^ 


>;(-.- 


--.r + E 


i^j 




Ikif^i) 


a 



[Xi - Va) 



(12) 



Using the properties of the roots of the Hermite polynomial 



j:in-r,)-\ E {n-r,)-' = L{L-l)/A, 



j(^i) l<i<j<L 

we thus conclude that the wavefunction (p is an eigenstate with eigenvalue 

TV 



E = -L{L - l)/4 + (1/2) E(A^<x - l)K. 



a=\ 



(13) 



(14) 



Although it is expected that this wavefunction is the lowest energy state in the subspace 
of fixed iVi, iV2, ■ ■ ■ , Nn, we were not able to prove it. However, in the case of SU{2), the 
small lattice diagonalization up to 8 sites confirms this conjecture. Moreover, the discussion 
below will also confirm this idea for the general case. For fixed number N^ of the electron 
number, the minimum of the energy is obtained when \N„ — iVg-'l is as small as possible for 
each pair o ^ o' . 

In the SU{2) case, the above result becomes 

E = (-l)L(L - l)/4 + {1/2)N^{N^ - 1) + {l/2)Ni{N^ - 1), (15) 

where N^ and N^ are the numbers of the up-spin electrons and the down-spin electrons 
respectively. For fixed number of electrons on the chain, i.e. for fixed N^, the minimum of 
the energy given in Eq. ([15|) is obtained when Sz = for even A'e, or when Sz = ±1/2 for 
odd A'^e. Therefore, the ground state is a spin singlet ( respectively spin 1/2 ) state for even 
( respectively odd) number of electrons on the chain. In particular, for an even number of 
electrons on the chain, the ground state energy is 

Eg = (-1/4)L(L - 1) + N^/A - NJ2, (16) 

while for an odd number of electrons it is 

Eg = (-1/4)L(L - 1) + i^^r. (17) 

The charge susceptibility of the ground state Xc is given by Xc^ = O'^Eg/ON^ = 1/2, 
independent of the electron concentration. Very unexpectedly, the charge susceptibihty 
is also finite at half-filling Ne = L, in spite of the existence of a metal-insulator phase 
transition at half-filling for this system. This is in contrast to the case of the periodic 1/r^ 
supersymmetric t-J model, where the charge susceptibility is divergent at half-filling, at 
which the metal-insulator phase transition occurs. 

To study the spectrum of the system away from half-filling, we follow the idea introduced 
in Ref. |l^. Let us define the operators 
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a] = -Kj +iqj, 

a. = Trj-iqj, (18) 

which satisfy the following commutation relations: 

[TTjjTTfc] = 

[qj,H] =iTij 

[K,,H] = -2t J2(q^-q,y\ (19) 

Then using the property of the roots of the Hermite polynomial we have 

[n„H] = -zq^, (20) 

and 

[aj,H] = —a] 

[a„H]=a,. (21) 

Therefore the operators Al{iy) = ajS^ ,i = 1,2, ■ ■ ■ , N^, where u = 0,±, z for the SU{2) 
case with Si = 1, will act as raising operators, while their hermitian conjugate Aj(z/) will 
act as lowering operators. It thus follows that the wavefunction 

Ne 

</>{"}>} = En(4(^.))"^^0 (22) 

p j=l 

with {n} = {rii, n2, ■ ■ ■ , nNj,ni > 0, {u} = (z/i, z/2, ■ ■ ■ , i^nJ, is either an eigenstate with 

energy 

E{n} = E + J2ni (23) 

or zero. 

Moreover, it is shown in the Appendix B that the operators J2i=i ^i{^i) with z/j = or z, 
and act annihilate the wavefunction 0, and also '}2i=i ^i{^)4'G = 0, confirming the conjecture 
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that is the lowest energy state in the subspace. We then arrive at the conjecture that the 
excitation spectrum of the system is of the form 

E{s) =E + s; 

sG (0,l,2,---,s^,,), (24) 

i. e., the spectrum of this t-J model consists of equal-spaced energy levels. Since the model is 
on a finite chain, Smax is finite. In the special case of SU{2), the small lattice diagonalization 
up to 8 sites suggests that the highest energy level is given by 

EmaAQ) = L{L - l)/4 - Q{Q - l)/2, (25) 

i.e., for an even number of electrons, 

Emacc = Eg + (l/4)iVe(4L - 3iVe), (26) 

where Q is the number of holes on the chain, N^. = L — Q is the number of electrons. For 
iVg = 1 or A^e = -^ ( half -filling ), this formular gives right results; moreover, at half- filling, 
this corresponds to all spins polarized in one direction. 

The feature of the t-J model spectrum consisting of equal-distant energy levels may also 
be seen by taking the strong interaction limit of the Sutherland-Calogero-Morse quantum 
system for a mixture of fermions and bosons, 

H = (-1/2) J: d'/dq^ + J: l'q^/2 + ^21(1- M,,)/{q, - q,)\ (27) 

i=l i=l i<j 

where there are N^. fermions with spins and Q spinless bosons, Mij permutes the positions 
of the particles i and j only. The mixture gas has equal-distant energy levels described in 
terms of effective harmonic oscillators. In the strong interaction limit, the elastic modes 
decouple from the internal degrees of freedom. Since elastic modes also consists of equal- 
distant energy levels, we thus are led to the conclusion that the spectrum of the internal 
dynamics, which is that of our t-J model, also consists of equal-distant energy levels. Further 
work is necessary for a fully complete proof that the t-J model full spectrum takes the form 
Eq. (H. 
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Finally, we would like to point out that the states of the t-J model in the whole Hilbert 
space are grouped into a structure of "spin supermultiplets" , as indicated by the small 
lattice diagonalization, similar to that of the periodic 1/r^ supersymmetric t-J model. Such 
pattern of the Hilbert space is related to the symmetries associated with the Hamiltonian. 
It is highly worth while to identify them more explicitly, and we would like to study these 
aspects in further work. 

In summary, a set of Jastrow eigenfunctions have been found for the t-J model, with the 
eigenenergies explicitly calculated. The expected ground state of the t-J model is included in 
this set of wavefunctions. The full spectrum of the t-J model is found to have equal-distant 
energy levels which are highly degenerate. It would be very interesting to understand the 
underlying symmetry principles that give rise to such simple Hilbert space structure. It 
remains to study various correlation functions, as well as the thermodynamics, for this 
strongly correlated electron system. It would also be very interesting to study the effective 
field theory for the low-lying excitations for this t-J model. 

We wish to thank Dr. James T. Liu and Dr. Nicolas Macris for conversations. In 
particular, we are very grateful to Dr. James T. Liu for his substantial numerical support. 
We also would like to thank the World Laboratory Foundation for the financial support. 



APPENDIX A: 

In this appendix, we provide a brief proof for the sum rule Eq. (^ for reader's con- 
venience. The same argument can also be found in previous works Refs. [|,1£]- Let 
X = (xi, X2,- ' ' , Xn)',Xi = X/xi, we wish to show 

1 



E-^ n 



i=l Xj&Xi -"^i -^i 



0,Vt = 0, 1,2, ■■■,(« -2). 



Indeed, the Vandermonde determinant V{X) has the property: 

where X = {xi,X2,- ■ ■ , Xn),Xx = {xi,X2,- ■ ■ , Xn, x). Therefore we obtain 

1 



E-i n 



i=l Xj&Xi -^^ ^J i 



Ex! ^^^^^ 



f, 'V{X,x,) 



I X 



\n—\ 



^^ ^ ^'ViX) ViX) ^ 



4 = 1 



Xn 



\x 



1 1 

Xi X2 

n-2 ri~2 



Xn 



1 

Xn 



xrv 



0. 



This thus proves the sum rule Eq. 



APPENDIX A: 



(Al) 



(A2) 



(A3) 



In this appendix, we shall show that the lowering operators Ai{z) and a^ give zero 
when acting on the Jastrow wavefunction given by Eq. (^). This will yield a partial 
confirmation that the wavefunction is the lowest energy state in the subspace where the 
number of particles of each spin component is fixed. For the b boson degrees of freedom, we 
have following property 



ai0 = O,iG {N, + l,N, + 2,---,L), 



(Al) 
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which is shown to be true using the sum rule Eq. (j^) and the property of the Hermite 
polynomial roots Z]j(^i)('"i~'"j)~^ = ''^j- The procedure to deal with the permutation operator 
Mij in Oj is very similar to that of proving (f) to be the eigenenergy state of the Hamiltonian, 
but we do not write the full details here. Combining the Eq. ( [AID with the fact Y^^=i ca = 0, 
we thus arrive at the following results: 

L 

^ a„0 = 0, 

a=N^+l 

E««</' = 0- (A2) 

Furthermore, we realize that 



A^e 



[J2A{z)]<P = 0. (A3) 

1=1 

We have been able to show this to be true, following the similar approach to handle the effect 

of the permutation operator Mij acting on the Jastrow wavefunction (p. In the particular 

case where N^ = Ni, the wavefunction function is a spin singlet and we may globally 

rotate the first identity of Eq. (|A3|) in the spin space, giving us 



Ne 



E(A(±)) 



i=l 



0. (A4) 



In summary, we have proved that it is impossible to construct non-vanishing eigenstates 
with the lowering operators and the wavefunctions in the subspace where the number of 
electrons of each flavor is fixed. 
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